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We present a self contained formalism modelled after the Brownian motion of a quantum harmonic 
oscillator for describing the performance of microscopic Brownian heat engines like Carnot, Stirling 
and Otto engines. Our theory, besides reproducing the standard thermodynamics results in the 
steady state enables permits us to study the role dissipation plays in determining the efficiency of 
Brownian heat engines under actual laboratory conditions. In particular, we analyse in detail the 
dynamics associated with decoupling a system in equilibrium with one bath and recoupling it to 
another bath and obtain exact analytical results which are shown to have significant ramifications 
on the efficiencies of engines involving such a step. We also develop a simple yet powerful technique 
for computing corrections to the steady state results arising from finite operation time and use it to 
arrive at the thermodynamic complementarity relations for various operating conditions and also to 
compute the efficiencies of the three engines cited above at maximum power. Some of the methods 
and techniques and exactly solvable models presented here are interesting in their own right and, in 
our opinion, would find useful applications in other contexts as well. 
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j3 I- INTRODUCTION 

In recentyears there has been an upsurge of interest in the interface between thermodynamics and quantum 
mechanics macro and nano/micro [2|], leading to a thorough reexamination of the basic concepts and principles 
of thermodynamics with ramifications in biological processes and soft condensed matter systems. New paradigms for 
notions of work, thermal machines etc. have emerged that provide deep insights into thermodynamics which in turn 
help enlarge its scope far beyond that envisaged originally and open up new possibilities Q • These developments force 
one to formulate questions concerning efficiencies of various heat engines using appropriate microscopic considerations 
Indeed, one has even started considering quantum heat engines which in principle have efficiency larger than 
Carnot efficiency though such cases require non equilibrium steady states. Such steady states can be reached via the 
use of the coherent laser fields or via quantum interference effects In a very recent experiment |3| realized a 

microscopic Brownian heat engine The most important ingredient both in the work of Blickle and Bechinger 

as well as in the proposals of Scully and collaborators [j| is the possibility that all the relevant parametrs can be 
very well controlled experimentally and thus the heat engine cycle can be precisely realized and it becomes desirable 
to have exactly soluble models of microscopic heat engines. 

In view of the way the experiments are carried out we need a fully dynamical model which should account for the way 
the system parameters like potentials or the external parameters like temperature are changed. Further the behavior 
of the engine should depend on various time scales foe example the time taken to reach equilibrium state. Such a time 
would depend on the scales of the damping in the system. Motivated by these requirements we develop in the present 
work an exactly soluble model of a microscopic Brownian heat engine. The model that we present is fully quantum 
mechanical. Our model enables us to examine many different possible experimental scenarios- (a) low temperature 
behavior where quantum effects are likely to dominate, (b) behavior under different relaxation conditions-for example 
the system could be underdamped or overdamped, (c) possibilities for the system to pass through nonequilibrium stages 
depending on the rate of change of the external parameters, (d) nonequilibrium conditions because the experimental 
time scales are smaller than the time it takes for the system to reach steady state. Our formulation is based on the 
Wigner function and quantum Langevin equations for an harmonic oscillator whose frequency is modulated in time. 
We also assume that the temperature of the environment is time dependent as well. These time dependences are 
needed to implement the heat engine cycle realistically. We calculate the time dependent Wigner function, all the 
fluctuation parameters and the entropy. These enable us to calculate thermodynamic quantities like work, heat and 
internal energy. 

A brief outline of the work is as follows. In Section II for later reference we briefly recapitualate the expressions for the 
efficiencies, both classical as well as quantum, for the three engines based on standard thermodynamic considerations 
with usual assumptions regarding the speed with which the various steps are carried out. In Section III we present a 
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qua ntum thermodynamic frame work based on Wigner phase space description for quantum systems which contains 
the classical framework as a limiting case. In Section IV we develop a general set up for computing various quantitities 
of interest and give two models of frequency modulation where the relevant equations are amenable to exact analytical 
results. In Section V we consider the situation when the diffusion constant is varied linearly and in Section VI analyse 
its ramifications on the efficiencies of Brownian heat engines. In Section VII we develop a systematic scheme for 
computing finite time corrections to the efficiencies of classical and quantum Brownian motors and then use these 
results in Section VII to examine the role they play in determining the efficiency of the Stirling engine at maximum 
power. Section IX contains our concluding remarks and further outlook. 

As noted above our working model for a heat engine is based on a quantum harmonic oscillator with a frequency uj 
interacting with a thermal bath at temperature T. Three typical engines that have been discussed extensively in the 
literature based on varying w and T appropriately are (a) The Stirling (b) Carnot and (c) Otto engines. Their cycles 
adapted to the harmonic oscillator model are schematically given below: 



Stirling Engine Q 



u>i,T c Isothermal W2,T C 
3 — > 4 

Tc 

Isochoric t I Isochoric 

Th 

2 i — 1 
ui\ , Th Isothermal u>2 , Th 



(1) 



Carnot Engine @,[TT| 



u! 2 > wi , T h > T c 



uj\,Th Isothermal W2,T^ 

1 — > 2 

Th 

Isentropic j" I Isentropic 

Tc 

4^—3 
u>4,T c Isothermal ^3,T C 



(2) 



Otto Engine [lz 



lo c ,T 2 Isentropic u h ,T h 
4 < — 3 

T\ 

Isochoric \. \ Isochoric 

T2 

1 — > 2 
uj c , T c Isentropic cuh , T\ 



(3) 



Lu h > uj c >, T h > T c , f3 c uj c = P\U) h , (3 h u h = fouc, ■ 

Here, the t's indicate the time taken to carry out the indicated step and (3 stands for 1/KbT. The calculations that 
we give in subsequent sections can be applied to any of these engines. 

The three prototype engines above thus involve suitable combinations of the following three steps: (a) isothermal 
i.e. w changing, T fixed or (b) isochoric i.e. u held fixed, T changing or (c) isentropic i.e both oj and T changing with 
ujjT fixed and one needs an appropriate formalism to compute the efficiencies under specific physical circumstances 
in which these steps are actually executed in an experiment. The present work has this as its major objective. 
Our principal results include (i) development of a self contained formalism for computing efficiencies of Brownian 
engines both in the classical as well as quantum conexts (ii) an exact analysis of the role of damping in the process 
of coupling the system to a bath at a higher temperature and its influence on the performance of the Stirling engine 
(iii) computation of the irreversible heat in isothermal processes and the derivation of complementarity relations (iv) 
a detailed analysis of the role of damping as well as finite time corrections on the efficiency of the Stirling engine at 
maximum power. 
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II. STEADY STATE EFFICIENCIES FROM THERMODYNAMICS 



To set the notation and for later reference we assemble here the standard thermodynamic considerations that 
enable us to compute the efficiencies for the three engines listed above both in classical and as well as in the quantum 
contexts. These are: 

1. the thermodynamic conservation law AU = AQ — AW; where AU: Change in the internal energy U ; AQ: Heat 
absorbed by the system; AW: Work done by the system 

2. AQ in an isentropic process a — > b = 0, 

3. work done in an isochoric process a — > b = 0, 

4. work done in an isothermal process a — > b = —[F(b) — F(aj] where F denotes the free energy of the system, 

5. the expressions for U and F for the harmonic oscillator: 

1//3, p~ { 



K B T 



u = 



fHo[n(co,T) + 1/2], n(u,T) 



(Classical) 
(Quantum) 



(4) 



F(w,T) = 



-ln(2sinh(/3?icj/2)) 



6. the expression for the entropy of a classical harmonic oscillator 
In an isothermal process one has 

AQ = AW = TAS. 
In the quantum case, for the thermal states p t h, 



(Classical) 
(Quantum) 



Pth 



(5) 



(6) 



(7) 



(8) 



Tr[e-^] 

where H denotes the hamiltonian for a quantum harmonic oscillator, one has for the von Neumann entropy 

S = K B [(n{u,T) + l)Ln(n(u>,T) + 1) - n{u,T) In n(u,T)]. (9) 

With this preparation we now proceed to compute the efficiencies of the three engines mentioned earlier both in the 
classical as well as quantum cases. These would then be compared with the results obtained from the microscopic 
theory developed later. 



A. Stirling Engine 

The efficiency rj g of the Stirling engine is defined as 

Work done by the system 



Vs 



Heat flow into the system at Th 



Classical 
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In the classical case the work done by the engine is given by 

AW^ 2 + AW3-4 = -[F(u u T h ) - F(uj 2 ,T h )} - [F(oj 2 ,T c ) - F(u u T c )] 

= K B (T h -T c )ln(^j, (10) 

and the heat absorbed at T h by 

AW^ 2 + AC/^ 2 + = -[F(w u T h ) - F(w 2 , T h )] +0+1(1--!-) 

2 Ph Pc 

= K B T h \n (^j + \K B {T h - T c ). (11) 
[ Note the factor of half in the second term on the RHS of the above equation. We will return to this later.] Hence 

^ cl = V - j-. Vc = i-^- (12) 

Quantum 

Proceeding as before and using the expressions for U and F appropriate to the quantum case, we have for the work 
done 

AWi_> 2 + AW 3 _>4 - -[F(u u T h ) - F(u2,T h )] - [F{u 2 , T c ) - F(u u T c )], 

. Kb t„ ■„ c^pm - k bT , h , w 

Vsmh(^fiwi/2)/ V slnn (Pc^i/2)/ 

and for the heat absorbed at T h 

AW^ 2 + AC/1^2 + AC/4^1 = -[^(wi.Tfc) - F(w 2 ,T fc )] 
+ [fiwi(n(wi,T h ) + 1/2) - hio 2 (n(cu 2 ,T h ) + 1/2)] 

+ 1 ([fiw2(n(w2,T h ) + 1/2)] - [^ 2 (n( W2 ,T c ) + 1/2)]) , (14) 



and hence 



ls 1 + Z/X ' 
, / sinh(^^ 2 /2) \ y= /^, ( sm\i{p c hw 2 /2) 
^wah^htiU!^))' /3 C ^ S inh(/3 c ?^w 1 /2) / , ' 



Z = 

2 



fiW2 



hui coth (fihtvjJx/2) ^{coth {(3 h hoj 2 /2) + coth (/3 c ftw2/2)} 



2 

In the limit /3u> << 1 rj^ goes over to the classical efficiency r)f as expected. 



B. Carnot Engine 

For the Carnot engine the efficiency defined as before 

Work done by the System 



Vc = 



Heat flow into the system at Th 
turns out to be the same in both classical and quantum cases and is given by 



5 



C. Otto Engine 

Here again the efficiency defined as 

Work done by the System 



Heat flow into the system during 2 — s- 3 
turns out to be the same in both quantum and classical cases and is given by 

to to U(3) - U(2) 

-('-£)■ « 17 » 

The expressions for efficiencies for the the three engines, realized here through a harmonic oscillator by appropriate 
changes of its frequency ( or equivalently its 'spring constant') and the ambient temperature, hold for idealized 
operating conditions as stipulated in equilibrium thermodynamics. These, for instance, demand that the isothermal 
changes of frquency involved in the Stiring or the Carnot cycles be carried out quasistatically i.e. so slowly that at 
each instance the oscillator remains in the state of equilibrium at that temperature and frequency. Such conditions are 
hardly ever met in practice and particularly in the light of the experimental work reported in [7[ there is an obvious 
need for developing a framework which brings into play aspects of approach to equilibrium , both in classical and 
quantum contexts, and is capable of furnishing a self-contained scheme for computing the efficiencies under realistic 
conditions. We develop such a scheme in the next section. 

III. EFFICIENCIES BEYOND THE STEADY STATE: A DYNAMICAL MODEL 

To go beyond the standard thermodynamic assumptions regarding the rate at which which various steps in a heat 
engine are carried out so that one can evaluate the performance of an engine under actual laboratory conditions we 
need a framework which treats the system modelling the engine as an open system and permitting proper inclusion 
of dissipative effects and the possibility of varying the ststem potential and the ambient temperature. In the present 
context, such a framework is provided by the dynamics of a quantum Brownian oscillator of frequency us in contact 
with a heat bath at temperature T.is described by the master equation [13[ 

|p = + \mu 2 q 2 ,p] - ^pW, T) + l/2)([g, [q, p]]) 

([$,{& P}]), (18) 

where q and p are denote the position and momentum operators obeying the commutation relations [q,p] — ih. 

For reasons given later it proves expedient to transcribe the quantum dynamics described by the master equation 
using the Wigner phase space description of quantum systems (3], [l5[ which associates with a density operator p a 
phase space function W(q,p) of classical variables q,p as follows: 



p^Wp(q,p) = Tr {pW(q,p)] 



W(q,p) = J dq' \ q + l -q'){q- l - q '\e l ^/\ (19) 

— oo 

We note here that we use the Wigner phase space description in preference to other phase space descriptions for 
two reasons: (a) it is the only one that that maps the quantum mechanical average of a product of two operators to 
the phase space average of the corresponding Wigner functions (b) its moments < q m p n > correspond to quantum 
averages of of the symmetrised operator (q m p n )s- For example < q 2 p > corresponds to the expectation value of the 
operator (q 2 p + qpq + pq 2 )/3. 

Use of the Wigner description turns the master equation into a Fokker-Planck equation for W(q, p) [l5| 



dq (m) dp ( dp 



d 2 

W(q,p,t), 



r,2 
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where 

1 



V(q,a) = -aq 2 , a = muj 2 , (20) 



and 



D = 2mhwK{n{oj, T) + -), n{w, T) = {e pHu - l)" 1 . (21) 

In the following the parameter a, the 'spring constant', will be taken to be controlled externally. 
The Langevin equations equivalent to the above FPE read: 

<?=-, (22) 
m 

p = -2Kp--^-V(q,a) + f(t), (23) 

< f(t)f(t')>=2D8(t-t'). (24) 
The Langevin equations (|2"2"|) — ([Mf lend themselves to a nice thermodynamics intepretation Q : Rewriting (|2"3")l as 

-(-2 K p + /(t))+p + |-F(g,a) = 0, (25) 



and multiplying it by dq and using 



one obtains 



<9V(g,a) 9F(g,a) 

dK = dg H da, (26) 

oa c*a 



- (-2kp + /(t))dg + d(p 2 /2m + F(g, a)) - ^'^ da. (27) 

oa 

The three terms in the above equation may now be identified in an intuitively plausible manner as: 

dQ = (-2np+ f(t))dq, dU = d(p 2 /2m + V), dW = - 9V ^ ,a ^ da, (28) 

oa 

leading to the energy balance equation: 

-dQ + dU + dW = 0, (29) 

with dQ ( -dQ) understood as the heat flow into of (out) the system and dW (-dW) as the work done by (on) 
the system. The stochastic averages of these quantities denote by dQ,dU and dW respectively relate directly to 
the corresponding thermodynamic quantities and capture the thermodynamic conservation laws. This self-contained 
approach is clearly more microscopic than thermodynamics as it provides a framework for computing not only the 
averages of these quantities but their probability distributions as well. 

We note here that while it is certainly possible to transcribe the master equation dynamics directly into equivalent 
quantum Langevin equations for the operators q and p but owing to their noncommutativity the crucial step (|2T[) 
needed to obtain a clear thermodynamic interpretation of the such Langevin equations would now involve terms like 
p.dq and would therefore be fraught with ordering ambiguities. 

The scheme described above for computing dQ, dU and dW together with the expression for von Neumann entropy 

S = K B [((r + l)ln(a+l)-alna], a = ^Det[V] - -, (30) 

for Gaussian states [16j ( which is what we would exclusively deal with) i.e. states p for which the Wigner distribution 
is a Gaussian: 



W(q,p) = — = exp 

V(2^) 2 Det[V] 



x Vx 
2Det[V] 



(31) 
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provide all that we need for the considerations below. Here V stands for the variance matrix, 

V=( <q2> <qP 2 > ), (32) 
\<qp> <p 2 > J' y ' 

and < • > denote averages with respect to the Wigner distribution. The uncertainty relations require that a be 
positive. Note that the set of Gaussian states contains the set of harmonic oscillator thermal states pth as a special 
case. 

Before proceeding further it is instructive to check that the Wigner description above together with the thermo- 
dynamic interpretation implied by (|28[) in the steady state does indeed reproduce the results given earlier for the 
efficiencies of the three engines using standard thermodynamic considerations. Thus, for instance, calculation of the 
efficiency of the Stirling engine involves computing AW i_>2 AC7 i_>2> Af/4_>i which in the present framework are given 
by 

AWi_> 2 = dW = muj <q 2 > T =T C duj, (33) 



UJ2 

/>4 fUJ2 

= / dW = - mcj < q 2 >T=T h dco, 

J 3 J wi 



(34) 



Further, from the FPE or the Langevin eqns it follows that in the steady state 

<p 2 >= D/2k, muj 2 < q 2 >= D/2mn. (37) 

These on using D = 2mKUj(n(uj.T) + 1/2) then give 



AWUa = - f 1 [»(w, T h + l/2]dw 



KBT h ln[ S ^ h ^2 )=F(l)-F(2), (38) 



v siiili(^wi/2)y 
AW 3 ^ 4 = - / [n(w, T c ) + l/2]dw 



= -^T c ln( ^^ )=[F(3)-F(4)], (39) 
Vsmh(/3 c wi/2)/ 

AC/x^ 2 = [«i(n(wi,T fc ) + i) - [ W2 (n( W2 ,T,0 + i)], (40) 

Atf 4 -i = i ^[waKwa.Tfc) + -)] - Mn(^ 2 ,T c ) + , (41) 

which are the same expressions as before and therefore one recovers the expression for efficiency given in Section II. 

( Note here that in computing AU^i we considered only the contribution from < p 2 > and not from < q 2 >, a 
question that will be examined in greater detail later.) 



IV. QUANTUM DYNAMICS UNDER TIME DEPENDENT CHANGES OF TEMPERATURE AND 

POTENTIAL 

We have seen in the previous section that the Langevin equations equivalent to the Fokkcr-Planck equation obeyed 
by the Wigner distribution lend themselves to a direct and transparent thermodynamic interpretation and that this 
interpretation in the steady state limit reproduces the standard thermodynamic results. To prepare ground for going 
beyond the steady state limit we now analyse the structure of the solutions of the Langevin equations at hand allowing 
for arbitrary time dependence in the potential and the diffusion coefficients and apply this framework to arrive at the 
exact solutions of the Langevin equations for three fairly realistic models. 
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The Langevin equations which in the present case are linear stochastic equations with additive noise may be solved 
to yield : 



fq(t) 
\p(t) 



= M(t) 



where 



9(0) 
p(0) 



M(t) 



solves the homogeneous equations 



From 



one has 



for the variance matrix 



d_ ( q(t) 
dt \ p(t) 



V{t) 



dt'M(t)M(t'Y 



u[t) v[t) 
mu{t) mv(t) 



1/m \ / q(t) 
-tow 2 —2k J [ p(t) 



< q 2 {t) > < q(t)(p(t) > 
<q{t)p{t)> <p 2 (t)> 







V(t) = M{t)[ V(0) + / dt'M-\t') 




2D(t') 



M (t ) ]M (t). 



(42) 



(43) 



(44) 



(45) 



(46) 



It is therefore clear that finding explicit solutions for the variances in situations where both ui and T depend on time 
depends on our ability to solve for M(t). We list below three physically meaningful cases where this is indeed possible. 

Case I: u independent of time 

For this familiar case the functions U(t) and V(t) which determine the matrix M (t) are explicitly given by 

(A+e- A -*- A_e- A +*) (e"*-' - e" A +*) 



«(*) = 



(A+-A-) 



-, v(t) = 



m(A+ — A_) 



X± = k ± \/ k 2 — i 



(47) 



Further, owing to time translation available in this case, we have M 1 (t) = M(—t), M(t)M(t') — M[t + t'), and 
simplifies to 



V(t) = M(t)V{0)M T (t) + / dt'M(t') 




2D(t-t') 



(48) 



Case II : uj 2 (t) = (l + §) , < t < T 

In this case investigated in fl7j | . the functions u(t) and v(t) in the range < t < T are given by 



u(t) 



/ + (t)/_(0)-/-(t)/ + (0) 
/ + (0)/_(0)-/_(0)/+(0) 
2 



, v(t) = m 



/+(t)/_(0) -/_(*)/+ (0) 

/V(o)/_(o)-/_(o)/+(o) 



f±(t) = e- Kt (t + a) 1 ' 2 J ±1/3 o &1/2 (^ + a ) 



X3/2 



1-— -, b 



T ' 



(49) 
(50) 



Case III : uj 2 (t) = u 2 e^ , < t < T 

In this case the functions / + (t) and /_ (t) in (|50|) are again given in terms of Bessel functions as 



f±(t)=e Kt J± a (ae^) , a 



2Tw 2Tk 
a = . 



(51) 



Having dealt with some exactly solvable cases where the frequency is changed in a specific way but the tempera- 
ture may be varied arbitrarily, we now illustrate how the formalism developed above lends itself to useful exact or 
approximate calculations leading to finite time corrections. 
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V. LINEAR VARIATION OF THE DIFFUSION CONSTANT 



We first consider the case in which the the harmonic oscillator with frequency uj is in equilibrium with a bath at 
temperature T characterized by a diffusion constant D . With uj held fixed, the the diffusion coefficient is changed 
linearly from its initial value D appropriate to temperature T to its final value Di appropriate to temperature T\ 
in a time r and then kept at that value thereafter. 



lu : held fixed 



D(t) 



£>i, t > T. 



(52) 



This situation pertains to the isochoric step in the Brownian engines and is relevant for discussions on aspects of 
decoupling the system from a heat bath at one temperature and recoupling it to another heat bath at a different 
temperature. 

For the case at hand, with V(0) chosen to be the variance matrix corresponding to the oscillator being at equilibrium 
with the bath at temperature appropriate to Dq 



V(0) 



/ A> 



2kui 






D 



2Km 2 Lu 2 



(53) 



we have from J48 



Using the relations 



Do ( u 2 {t) 



■v\t) 



2k \m 2 u> 2 

<q(t)p(t)>=2 [ dt'v(t-t')v(t -t')D(t') 



dt'v 2 (t-t')D(t'), 



<p 2 (t) >= 



mD ( u(t)u(t) 
2k 



+ v(t)v(t) \ +2m 2 dt'v 2 (t-t')D(t'). 



(54) 
(55) 
(56) 



1 1 d 



2K2dt V m 2 uj 2 



1 d 



2k dt \ m 2 ui 2 



1 1 d 



vv 



2K2dt \ m 2 uj 2 



(57) 
(58) 
(59) 



which follow from 



we obtain for t > r 



-2kv H ; u 

m 



-mu) 2 v, 



(60) 



< q 2 (t) >= a(t) < q 2 > +(1 - a{t)) < q 2 >j 

< p 2 (t) >= p(t) < p 2 >o +(1 - P(t)) < P 2 >l 



(61) 
(62) 



where 



a(t) = 

m) - 



dt'[m 2 uj 2 v 2 {t')+u 2 {t% 



t-T 
t 



dt'[m 2 v 2 (t') + 



u 2 (t'). 



(63) 
(64) 
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In the limit t — >• oo limit both a(t) and go to zero and hence < q 2 > and < p 2 > assume their respective 
equilibrium values. The parameters a(t) and /3(t) thus interpolate between the initial and the final equilibrium values 
of < q 2 > and < p 2 > and quantify the approach to equilibrium. In the following we consider the case when t = r 
i.e. the situation that obtains immediately after the bath has reached the state characterized by the final value of the 
diffusion coefficient. Evidently as far as the system is concerned we are dealing here with a non equilibrium state as 
the system has not yet had time to equilibriate with the 'final' bath. 

Putting t = t in (163[) and (|64[) and denoting a(r) and /3(r) simply as a and j3 we obtain on substituting for u and 
v from (|T7| and carrying out the relevant integrals 



1 



13 = 



(x~y) 2 
1 

(x- yf 



(x + y) 1(1- 



-2y 



e -2«))_J^(l_e-W; 

/ (z + y) 

4 ^ n _ e - (»+!/) 
(x + y) 



(65) 
(66) 



where a; = [k + \J n 2 — uj 2 ]t, y — [k — \J n 2 — w 2 ]r. We now examine the behaviour of a and /3 in the overdamped and 
weak dissipation regimes respectively. 



Overdamped Case 

2 

In the overdamped regime i.e./t >> w, x sa 2kt >> 1, y s» << 1 one finds that 

1 - e~ 2 « 



22/ 
1 - e- 2 * 
2x 



1 as y — » 0, 
for a; >> 1. 



(67) 
(68) 



Underdamped Case 

On the other hand, in the weak dissipation regime k << uj,x w k + ico,y m k — iuj and we have in the limit kt — > 

'(l_ e -2«r) 



2kt 

1 - KT | 1 



sin cjt 

2\ 



(69) 
(70) 



(1 — e 2KT ) / Sill a'i 

- KT 



2kt 



1-KT 1 



smwT 



CUT 
2\ 



(71) 
(72) 



Note that a > (3 in both the cases. In fact this is always true - it can easily be shown that with x, y defined as before 



a-0 



x + y 



dt[e "" -( 



2xt _ p -2yf\ 



x - y Jo 

and hence a > /3 by virtue of the fact that the integrand is always positive. 



(73) 



VI. EFFECT OF TIME SCALES ON EFFICIENCIES OF BROWNIAN MOTORS 



We recall that in the calculation of the efficiency of the Stirling engine from standard thermodynamic considerations 
presented in Section II, we had drawn attention to the factor of 1/2 in the expression for AC/^i. Likewise in the 
computation of the Stirling engine using the quantum stochastic thermodynamics in the steady state we had noted 
that in computing only < p 2 > /2m contributes to and not muj 2 < q 2 > /2. This seemingly ad 
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hoc prescription can now be understood at a deeper level in the light of the analysis in Section V leading to eqns. 
(foTj) - (|64[) . It is clear from the discussion therin that, in general, the expression for AU^i should be taken to be 



ACW = (1 - 0) 



< p 2 >i < p 2 > 4 

2m 2m 



+ (l-a)mw I J, (74) 

where a and /3 depend on various time scales involved. Indeed in the overdamped regime a — > 1 and (3 — > and one 
recovers the earlier results. The mystery behind the factor of 1/2 in (ITT1) and that behind retaining the contribution 
from < p 2 > /2m in ([36]) alone is thus resolved. On the other hand, in the weak dissipation regime where both a and 
(3 are close to 1 the situation is very different and this has significant consequences for the the relative magnitude of 
classical and quantum efficiencies under same operating conditions as discussed later. Further, since for a harmonic 
oscillator < p 2 > /2m = muS 2 < q 2 > /2 we can rewrite the above equation as 

AU ^ 1=2 »{^2m 2m-) > " = 1 " — ' (75) 

Using this expression in the calculation of the classical and quantum efficiencies for the Stirlng engine given earlier 
respectively we obtain 

^ cl = Vc f x > (76) 

1 + rjcfi In - 

and 



i-yy* 
Vs i + z/x ' 

/ sinh(/3^ 2 /2) \ Y= fh, / sinh(/3 e ^ 2 /2) \ 
n ^sinh(^fi Wl /2)y' ' (3 C I1 ^sinh(/3 c ^ 1 /2)y' ' l " J 

Z = y [fiwi coth (fi h hu)i/2) - hjj 2 {(l - (J,) coth ((3 h hu 2 /2) + ;U coth( (3 c hu 2 /2 )}] . 

The appearance of the parameter /i appearing here may be viewed as a phenomenological way of incorporating non 
equilibrium effects arising from decoupling of the system from one bath and recoupling it to another. 
In terms of dimensionless quantities a, b, c as 

/3 c fiwi = a, — = b, ^ = c (78) 
the expression above for the efficienies 77^ and r]f in the classical and quantum read 

Vci(P,c) = 1 ~ C (79) 

1 + M- 



ln6 



r? g (a,6, c) 



V sinh(ac/2) / V sinh(a/2) J 



/sinh(a6c/2)\ ac , . , . a&c,, „ , , , , . 

ln -^—T? — + TT coth(ac/2) — ((1 - u) coth(a6c/2) - /icoth(a&/2)) 

smh(ac/2) / 2 2 



(80) 



In the experiments of Blickle and Bechinger Q 

a = 9.50065 x 1CT 7 , 6 = 2.04922, c = 0.845272. (81) 
With b fixed at this values we plot below the ratio R — r]%/r)f as a function of a, c for two representative values of /i. 




Fig.l The ratio R = rf^/rf^ plotted as a function of a = j3 c TvjJ\ and c = j3h/ C for 1^2/^1 = 2.04 
and (a) fi = 0.01 (b) fi = 0.5 

VII. FINITE TIME CORRECTIONS: COMPLEMENTARITY RELATIONS 

We next consider the situation when the system starts out at equilibrium with a bath at temperature T, and the 
frequency is changed from its initial value ujq to its final value u)\ in a finite time either isothermally (T held fixed) or 
isentropically (lo/T held fixed) and focus on computing finite time corrections to the standard thermodynamic results. 
Referring to (p])-© we recall that while Stirling and Carnot engines involve the former operation, Carnot and Otto 
engines involve the latter. The scheme for computing finite time corrections developed below is similar in spirit to 
the adiabatic approximation in quantum mechanics and is a variant of the method formulated in [18[ in the context 
of the Fokker-Planck equation adapted to the equations for the moments themselves. 

The equations for the second moments that follow from the Langevin or the Fokker-Planck equation may be written 

as 

±X(t)=A(t)X(t)+Y(t), (82) 



df 



where 




<q 2 >\ I 

X(t)= \ <qp> \ ,A(t)= [ -mu; 2 (t) -2k ± \ ,Y(t) = \ ). (83) 
<p 2 > J \ -Imu? -4k 

(At this stage, as indicated, we allow the frequency and the diffusion coefficients to be independent functions of t, 
Later however, we would specialise to situations appropriate to isothermal or isentropic variation of the frequency.) 
Putting t = st and expanding X(t) as 

X(t)=X^(s) + -X^{s) + -- - , (84) 

T 

we obtain 

A(s)X (0 \s) + Y(s) = => X^(s) = -A- 1 (s)Y(s), (85) 

XW( S ) = A- 1 ( S )^-X^( S ). (86) 
as 

The first of these equations can be taken to describe the situation where the system is in the steady state corresponding 
to the instantaneous values of ui and D and the second as describing deviations from this steady state. These equations 
then give 

< <z 2 M > (0) = 2 J^ {sW < «(-)?(«) > (0) - o; < p 2 (s) > (0) = (87) 



1 
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and 



<q 2 (*)> 



8k 2 + 2uj 2 (s) d 



8kuj 2 (s) ds 
1 d 



muj 2 {s) ds 



< q 2 {s) >(°) 
< q(s)p(s) >W + 



4Km 2 w 2 (s) ds 



to d 
~2~ds 

,2, ,2 



< q 2 (s) > 



(0) 



00 > (1) = -I 1 ^^^ < <? 2 (*) > (0) < P 2 ( S ) > (0) ] 



4k 



(89) 
(90) 



These equations together with (|87|) give finite time corrections to the variances. As the diffusion coefficient is a function 
of both uj and T we now specialize to the situations where (a) uj is time dependent, T is held fixed (Isothermal Case) 
(b) oj,T both are time dependent but ui/T is held fixed (Isentropic case). With this in mind we may rewrite expression 
for < q 2 (s) >^ which we would need shortly as 



< q 2 (s) >« 

h 

Annuo 2 



4k 2 \ ( , . 1\ doJ {An 2 \ d ( . . 1 



ds 



ds 



(91) 



In the isothermal case both the terms on the RHS contribute. On the other hand in the isentropic case only the 
first term contributes as during this process uj/T and hence n(w,T) are held constant. In the following we confine 
ourselves to the isothermal case and give the results for two physically interesting limiting cases corresponding to the 
overdamped and weak dissipation dissipation regimes. 

1. Overdamped regime( k >> uj) 

< q 2 (s) >« = 

'0h 



muj 5 (3 



2. Weak dissipation( k << uj) 



uj\ fBhuj\ f0fvjj\ 

-M-) +(-) 



cosech' 



/f3huj 



{to)' 



(92) 



2 ^ ^(i) = 



< <Z» > 



{~ 



cosech | — — | 

v 2 ; 



I 2 



duj\ 
~to) 



2ktocj 3 /3 

We now compute expressions for the irreversible heat Q m : 

Q ixi = F(f)-F{f)-/±W i ->f, 
in an isothermal process at temperature T from i — > f arising from finite time corrections. Recalling that 



(93) 



(94) 



AW^ f 



muj < q > duj 



and that 



we have 



< q 2 >(°) duj = F(i) - F(f), 



1 f 10 *' Yj 

Q iTT = - mui < q 2 > (1) duj = T— 

T ./„,.. T 



(95) 



(96) 



(97) 
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From the way the quantity £ is defined above it is clear that it would, in general, depend on both T as well as on the 
manner in which u) is varied from its initial value Wj to its final value ujf in the time r. We now turn to the question 
as to what would be the minimum value of Q; rr and hence that of £ in the weak dissipation and overdamped regimes. 
Using (p?2"j) and (j9"3")l and parametrizing w(s) such that w(0) = u>i, to(l) — ujf we have 



1? Jo dag(u(a)) ( ^ ) 



fir 



duj\ ' 



ds ) 



L- £ dsg(u(s)) ) 



2k(3t 



\ds J 



(K » id) 
(k « LU). 



(98) 



where 



1 

2^* 



f/3huj\ , f /3hoj\ f0fvjj\ 
' 1 coth 1 11 ! 



V 2 



cosech 



{—) + [-) 

fihw 



cosech 2 



(k » lo) 



(99) 



The expression Q m involve the functional 



I[uj] = 



(K « UJ) 



did 



dsg(ui(s)) ( — 



If we define fl = /(w) and its inverse as uj = h(Jl) and choose f(uj) to satisfy 



4f(w 



du) 

then we find that curve uj(s); < s < 1 which minimizes /[w] is given by 

w(*)=M/(w(0))(l -*) + /(("(!))*)), 

and its minimum value of by 

/mi„H = (/(^(l))-/(^(0))) 2 . 

In the classical limit (f3huj <<< 1) in the strong damping regime (k >> lu) we have 

ff(w) = ^=>/(w) = -i; /*(n) = -i, 



and hence 



min _ 2nK B T 
Win ~~ 



w(l) w(0) 



(100) 

(101) 

(102) 
(103) 

(104) 
(105) 



in agreement with the results of Sekimoto and Sasa [18| ( modulo an erroneous factor of 1/4 in the value of Qf^as 
quoted). This minimum value is realised along the curve 



w(0) w(l) 



;0 < s < 1. 



(106) 



Thus in the classical limit, in the overdamped regime, we obtain the following inequality for the product of the 
irreversible heat and the time taken to execute the step : 



t x Q irr > 2kK b T 



1 



1 



cj(1) w(0) 



(107) 



Such a relation is referred to in the literature as a thermodynamic complementarity relation, an analogue, both in 
spirit and form, of the energy-time uncertainty relations in quantum mechanics. It should however be noted that the 
RHS of the above complementarity relation is independent of K. 
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Again in the overdamped regime, but now in the low-temperature limit i.e /3fijj —> 0, we find that g(u) ~ /3h/Aui 3 
and the complementarity relation becomes 



r x Q irr > H(2k) 



1 



1 



-1 2 



vMl) vMo) 



(108) 



and one finds that hbar now does appear on the RHS as one would expect in the limit of low temperature where 
quantum effects become significant. 

Turning to the weak dissipation case, one finds that expression for g(u) is such that the relevant integral in (|101[) 
can be given in a closed form leading to the following complementarity relation: 



T X Qi rr ^ 



K B T 
2k 



Aarm(/3fcj(l)/4) \ 
° 8 ^tanh(/3/ia;(0)/4) J 



valid for all values of T. In particular, in the classical limit it becomes : 



t x Q irr > 



K B T 
2k 



log 



"(0) 
W (l) 



and the curve ui(s), < s < 1 which minimizes Qin- now turns out to be 

w(s) =o;(0) (1 - s)/2 o;(l) s / 2 . 



(109) 



(110) 



(111) 



We emphasise here that though we have presented explicit expressions for Q m for the weak and stong damping 
regimes, the results in (|87|) and (|88| enable us to derive expressions for Q lrr both for classical as well as quantum case 
without any specific assumptions on relative magintudes of k and uj. Further, in the classical limit we find that for 
an isothermal process from i — ¥ f carried out in a finite time r, Q m has the structure 



hrr = T-, 
T 



(112) 



where E is independent of T. On using the fact that for an isothermal process F(i) — F(f) = TAS^f, we may 
rewrite (l94l) as 



AWi 



T AS.. 



(113) 



This provides a convenient and physically useful way of parametrizing deviations from quasi-staticity in that in the 
limit Ti^t f —¥ oo one recovers the familiar results of equilibrium thermodynamics. 

VIII. EFFICIENCY OF THE STIRLING ENGINE AT MAXIMUM POWER 



In this section we would use the results of the previous section to analyse the efficiency of the Stirling engine at 
maximum power very much in the spirit of the earlier works in the context the Carnot cycle. We would closely follow 
the works of Schmiedl et al @ and of Esposito et al [ll[ who analysed the question of the efficiency of the Carnot 
cycle at maximum power in the limit of low dissipation from fairly general considerations. In particular, in [TTj it 
was shown that the Carnot efficiency at maximum power 77* is bounded below by rj c /2 and above by r\ c j (2 — rj c ) and 
that while the Curzon-Ahlborn efficiency (l9| is reached in the limit of 'symmetric' dissipation, the upper bound is 
realized in a completely asymmetric limit and coincides with the universal upper bound derived in [20j and [2l| from 
somewhat different considerations. 

Consider the situation when the isothermal steps 1 — > 2 and 3 — > 4 are carried out in finite times and r c 
respectively as indicated in ([TJ. Power generated during the Stirling cycle is then 



P 



AWi- 



AW 3 



Th 



(114) 



Also, as we have seen in Section VI, that the expression for the efficiency for the Stirling engine can be written as 

AWx^i + AW 3 ^4 



fiK B (T h - T c ) + AWl 



(115) 
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where /j, ~ in the weak coupling regime and equal to 1/2 in the overdamped regime. 

Using (TITS)) and putting E^ 2 = E 3 ^ 4 = E c , ASi_> 2 = -AS 3 ^ 4 = AS (fTT4"j) and (fTT5)) become 

p _ (Tfe - T C )AS - T h Y, h /T h - T c Y, c /t c 

Tc + T h 

d _ (^j - T C )AS - Th^h/rh - T c Tj c /t c , , 

7/s ~ pKBiTh-TJ+ThAS-ThXh/Th- 1 j 

Maximizing P with respect to t/j and r c one finds that P attains its maximum value for 



Substituting these values for th and r c in (| 1 1 7[) one finds that the efficiency for the Stirling engine at maximum power 
is given by 



* i 1 + \[W) 



nf* = ^ v * n ~* J . (120) 

(i + \/*&) + ft(i-fc) + i% 

We now consider two cases: 
Case A /i = 

In the extreme weak dissipation regime i.e. [i — 0, one recovers results similar to those in @,[lH in the context of the 
Carnot cycle : 



1. In the symmetric case i.e E c /Efr = 1, vf* equals the Curzon-Ahlborn efficiency t\ca — 1 — \/T c /Th'. 

^ = 1 : vt* = VGA- (121) 

2. -qf* is bounded by rj c /2 and r) c /(2 — r/ c ) 

Vc/2 < vf* < Vc/(2 - Vc). (122) 
The upper and the lower bounds respectively correspond to E c /E^ — > and E c /E^ — > oo 

Case B fi ^ 

For small but non zero [i <\ log(o;2/w 1 ) these results get modified to those given below 

1. In the symmetric case i.e E c /E/, = 1, i]f* is less than the Curzon-Ahlborn efficiency [lj| t\ca = 1 — \jT c jTh: 

, c l* 

^r- = 1 : '7 S = 7 V < ( 123 ) 

E/j 1 I ( M W 2t)c A \ 

^iog(cj 2 /"i)y 

2. j^ 1 * is bounded by r? c /2 and rj s /(2 — rj s ) 

Vc/2 < r)f* < 77,-/(2 - r h ). (124) 
As before, the upper and the lower bounds respectively correspond to E c /E^ — » and E c /Eh — > oo 
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On the other hand if fi > \ log(oj2/wi), one finds that 

rjf* < Vc /2. (125) 

In the figures below we display the bounds on r]f* for \i = 0.001, 0.1, 0.2, 0.4 with ioijw\ taken to be 2.05 where we 
also give the plots for r\ c ,r\cA and rj c /2 for comparison. 




Fig. 2 Efficiency r)f* of the Stirling engine at maximum power as a function of the 
Carnot efficiency r] c . The graph 77° versus rj c , for all values of E c /E/i, lies in the 
shaded regions given here for U2/0J1 = 2.05 and (a) \x = 0.001 (b) \x = 0.1 (c) \i = 
0.40. 0.2 and (d) \x = 0.4. The graphs of rj c , rjc a and rj c /2 versus r/ c are displayed 
for comparison. While in (a)-(c), fi < h log(u;2/u;i), (d) corresponds to the case 
when ix > i log^/wi). 



IX. CONCLUSIONS 



In this work we have developed a microscopic framework for computing efficiencies of quantum/classical brownian 
motors realized by a harmonic oscillator. Two exactly solvable models for frequency modulation are presented. In 
the context of the Stirling Engine we have analytically treated the question of coupling the system at equilibrium 
with a bath at low temperature to a bath at higher temperature and the role that various time scales play in this 
process and have shown that these have strong influence on its efficiency. Further we have developed a procedure for 
computing finite time corrections to the quantitities of interest needed for calculating the efficiency of the the three 
engines considered - Stirling, Carnot and Otto, and have derived the thermodynamic complementarity relations in 
the overdamped and underdamped situations both in the high as well as low temperature limits. In the spirit of the 
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works reported in Q and [TT| on the Carnot engine, we have analysed in detail the efficiency of the Stirling engine 
at maximum power and have investigated the role of dissipation parameters thereon. Though in this work we have 
exclusively considered interactions between the system and a thermal bath, the formalism can easily be extended 
to situations where the thermal bath is replaced by a squeezed thermal bath bringing with it new parameters and 
thereby ushering in new possibilities that have no classical analogues. 

In the present work we have modelled the three heat engines after the quantum harmonic oscillator. It is of 
interest to carry out a similar analysis for finite state quantum systems. Indeed the entire perspective on heat pumps, 
refrigeratators and heat engines developed in [3] is based on the simplest of quantum systems - a qubit. Though in 
that context a convenient Fokker-Planck framework would no longer be available we expect that the methodology 
developed here applied directly to the master equation would be useful there as well. We hope to return to this and 
related questions in the near future. 
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